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, A procedure to derive bounds on coupling strengths of exotic particles to nucleons from the 

neutrino signal of supernovae is outlined. The analysis is based on a model independent 
calculation for the emissivities for the exotic, detailed simulation for the evolution of the early 
proto-neutron star as well as a Likelihood analysis. As an example we derive confidence levels 
for the upper bound of the size of gravity only extra dimensions. 



In the aftermath of a core collapse supernova we find a system described by extreme param- 
eters: densities of several times nuclear matter density with simultaneous appearance of tem- 
2 ■ peratures of several tens of MeV. Indeed, the matter is that dense and hot that even particles 

, interacting as weakly as neutrinos — the by far most efficient conventional cooling mechanism — 

get trapped. Therefore, the cooling of the nascent proto-neutron star (PNS) has to happen 
on diffusion times scales of tens of seconds. It is this delayed cooling in the standard scenario 
^ , that makes core collapse supernovae one of the most efficient elementary particle physics labs 

possible: any additional cooling mechanism will influence the time structure of the neutrino 
signal as long as its interaction with the medium is not too strong that it gets trapped as wellEl. 

To make this kind of an argument more quantitative several ingredients are in principle 
necessary: first one needs a reliable method to calculate the primary emission rates of the 
exotic under discussion in the inner core. The calculated emissivities are then to be fed into a 
simulation capable of describing the first several tens of seconds of a nascent PNS. The results of 
these simulations carried out for different coupling strengths of the additional cooling mechanism 
are the to be compared with the data and analyzed with the appropriate statistical method. This 
leads then to bounds on the coupling strength of the exotic studied. The chain of arguments is 
illustrated in figure (|T|). Also indicated in the diagram there is a possible short-cut, here labeled 
as Raffelt criterium. This particular criteriumlll is one of a group of criteria meant to allow to 
read a constraint off the emissivities directly. Although useful to get a very rough estimate of 
the allowed strength of an exotic one should keep in mind that the emission rates of particles 
in the nuclear medium are strongly temperature and energy dependent. The core temperature 
in the supernova, however, depends on whether or not there is an additional cooling mechanism 
present. Thus it is difficult if not impossible to fix a fudicial temperature and density to be used 
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Figure 1: The structure of the analysis. Between input (framed in rounded boxes), calculated quantities (boxes) 
and analysis tools (ellipses) is explicitly distinguished. The right column shows the comparison of the results of 
the simulations to the measured integrated number of counts (solid: a2=0, short dashed: a2=0.01, long dashed 
02=0.1) as well as different contours of the log-Likelihood function q. 



in the criteria that is vahd for all kinds of particles (note that the temperature dependence of a 
particular production process depends on the energy /momentum dependence of the fundamental 
coupling vertex of the exotic to the medium and thus depends of the properties of the particle 
studied). Thus, to get a reliable bound, a detailed simulation is required after all. 

Let us now go through the structure of the analysis in more detail: the first step is the 
calculation of the emissivity of the exotic. In our work we follow the method developed in ref. i. 
It is based on the observation that, if the radiation is soft, the production rate can be expressed 
in terms of the on-shell nucleon nucleon scattering data directly — the process is dominated by 
radiation off external legs. Obviously this kinematic requirement is best fulfilled in degenerate 
nuclear matter as it occurs in the late cooling phase of a PNS. It was shown in ref. Bl that our 
model independent approach leads to rates for neutrino-anti-neutrino bremsstrahlung that are 
significantly smaller than those used so far. This finding was confirmed recently, as reported 
in ref. i. The separation of scales is not that clear when we move to the physical conditions 
present in the first 20 s after the collapse. With temperatures of some tens of MeV and densities 
of several times nuclear matter density, matter should he expected to be neither degenerate 
nor non-degenerate. However, the investigations of ref. a suggest, that the relevant expansion 
parameter, namely the radiation energy in units of the average nucleon-nucleon energy, still 
is smaller than one even here. Using this approach emissivities where already calculated for 
axions Bl as well as gravition only extra dimensions (GODs) In addition the latter results 
where used as a basis for bounds derived on the ground of cosmological arguments 10. 

Independent of the production mechanism, it is always possible to cast the production rate 
in the following form 

MeV/ba,Ws^ (1) 

To be concrete, let us concentrate on the production of gravitons in the presence of n gravity 
only extra dimensions (GODs) 0. There we find p2 = 5.42 , 02 = 5.1 x 10^i?2) = 6.5, 



03 = 1.4 X 10^^i?3. Here, ns is the baryon number density, T is the temperature. The size of 
the extra dimensions, Rn, is given in mm. The emissivitv depends weakly on the composition 
of matter through the function xi^n^ Xp), defined in ref. Q. 

Once the emissivities are derived we have to quantify their influence on the neutrino signal 
of SN1987A. The by far cleanest way to do so is to implement the emission rate of eq. (|^) in 
a simulation of the early PNS. For the present study we use the codes described in ref. B. The 
equilibrium diffusion approximation employed in the numerical code provides a fair description 
of the total neutrino luminosity and its time structure (c.f. figure (|^)) and thus forms a sound 
basis for the present investigation. However, in order to deduce the anti-neutrino spectrum an 
additional assumption is needed, namely that the neutrino emission happens from a common 
neutrino sphere whose temperature is determined by the requirement for the optical depth to 
take some fixed value. In order to minimize the influence of this assumption on the bound 
derived we studied a variety of values for the neutrino sphere temperature around the central 
value. In addition, the mass of the PNS very strongly influences the neutrino luminosity t3. 
Since we only know, that the baryonic mass is bounded in the interval IAMq < M < 2.OM0, 
we studied the whole mass range. 

Having neutrino emission rates at hand we now have to compare to the data and develop a 
ol that quantifies the agreement. Thus we need to define a Likelihood function. Following ref. 



we use 
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where N£) is the total number of observed neutrino arrivals, and tf, tf*! • • • the times at which 
neutrinos actually arrived in the detector D. Nu^t) is the total number of neutrinos that the 
model with assumptions /, PNS mass M, anti-neutrino temperature Tp^, and exotic-coupling 
Un predicts will have arrived in the detector D up until the time t. At can be any interval small 
enough that the probability of detecting more than one count in any one bin can be taken to be 
negligible. Ultimately it will be absorbed into an overall normalization constant. 

The likelihood function provides us a with quantitative tool to compare models of PNS 
evolution. Since we treat not only the coupling to extra dimensions, a„, but also M and Tp^, 
as parameters, the likelihood function is a function in a three-dimensional space. This function 
has a minimum at a baryonic mass of M = I.SMq, 02 = and Tp^ = 1.1T°. T° denotes the 
"reference anti-neutrino temperature" defined above using the optical-depth prescription. We 
can now assess all other models by looking at the log of the ratio of their likelihood of a particular 
model to this "most likely" model: 

I ^ N 1 C{{data}\an,M,Tp^,I) \ 

q{a^, M, r,J = - log (^^(|rf,,„}|o,i.5Mo,l.iro,/)j 

The function q is then a function in this same three-dimensional space. 

For the case of two GODs we discuss first the situation where the extra dimensions have zero 
radius, and so 02 = 0. Contours of q in the resulting two-dimensional x — M plane are shown 
in the left panel of the contour plots in Fig. (|l|). Here the horizontal axis is the temperature 
normalized to the reference anti-neutrino temperature, x = Tp^/T°, which is varied to explore 
the sensitivity of the results to deviations of the neutrino temperature from that obtained using 
the optical depth prescription. We see that although M = l.SM©, T = 1.1T° is indeed a 
minimum of the function g it is a rather shallow minimum: varying the neutron-star mass and 
the anti-neutrino temperature over the range considered here does not reduce the likelihood 
greatly. Ultimately this refiects the weakness of the constraint that the SN 1987a data provides 
for these parameters. 



The situation is rather different as we move away from 02 = 0. In the middle panel we 
display similar likelihood contours in the M — log(a2) plane for the case Tp^ = T°, while the 
right panel shows contours in the x — log(a2) plane for the baryonic mass M = I.GMq. These 
panels show that the likelihood function decreases rapidly for 02 > 10~^ — regardless of the values 
of the poorly-known parameters M and Tp^. Such large values of 02 are essentially two orders 
of magnitude less likely than the "most likely" model. This is a contrast to smaller values of 
02, where the differences in likelihood are comparable to those seen in the M — Tp^ plane. Thus 
statements about the most likely value of 02 in this smaller-a2 regime will be sensitive to the 
ill-constrained information on these neutron-star parameters. Consequently we will not make 
any such statements here. However, the two rightmost panels give us confidence in our ability 
to derive a bound on 02 — as opposed to a most-likely value — since it is clear that certain values 
of this coupling can be well-excluded, completely independent of details of the PNS modeling. 
To actually derive a bound we now integrate (or "marginalize"! over all possible values of M 
and Tp^, using appropriate weight functions (for details c.f. ref. S). We find that the possibility 
that there are two compact extra dimensions with radii larger than 0.66 /xm is excluded at the 
95% confidence level — as is the possibility that there are three compact extra dimensions larger 
than 0.8 nmi. 

To summarize, in this talk a scheme to deduce bounds with a well defined statistical meaning 
for the coupling strength of exotic particles to nucleons was presented. The method is quite 
general and it was demonstrated that the bounds derived do not depend on parameters poorly 
constrained. However, improvements are possible, and should be studied. For instance, many- 
body effects are expected to modify the emissivity. We are working on this problem. 
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